Entanglement trapping in a non-stationary structured reservoir 
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We study a single two-level atom interacting with a reservoir of modes defined by a reservoir 
structure function with a frequency gap. Using the pseudomodes technique, we derive the main 
features of a trapping state formed in the weak coupling regime. Utilising different entanglement 
measures we show that strong correlations and entanglement between the atom and the modes are 
in existence when this state is formed. Furthermore, an unexpected feature for the reservoir is 
revealed. In the long time limit and for weak coupling the reservoir spectrum is not constant in 
time. 
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I. INTRODUCTION 

In recent years, entanglement and quantum correla- 
tions have attracted the attention of many physicists 
working in the area of quantum mechanics [l|, [2| . This is 
due to the ongoing research in the area of quantum in- 
formation [3, and also because of the advances made in 
different experimental disciplines .such as in ion traps [4( 
and Bose-Einstein condensation [5|, |6[ . Developments in 
the field of cavity QED, where experiments in the strong 
coupling regime are carried out 0, H[ provide plenty of 
motivation for studying quantum information and en- 
tanglement. Theoretical studies are also important in 
the context of atom-light interactions inside structured 
reservoirs [9] such as resonant cavities or photonic band 
gap materials. The theoretically predicted atom-photon 
bound state could also lead to entanglement and this can 
also be linked to another problem: that of atom-laser out- 
coupling from Bose-Einstein condensates [Iol-[l2j, where 
analogous effects were predicted in the past. 

When quantifying entanglement between an atom and 
a reservoir of modes, the modes can be treated collec- 
tively [13| . The system is described in terms of two sub- 
systems and one can use existing bipartite entanglement 
measures. This of course does not permit the study of 
entanglement between individual reservoir modes. It is 
also possible to partition the reservoir and then quan- 
tify entanglement between different parts of the reservoir 
[14{ . A different approach is that offered by a recently 
proposed measure, the density of entanglement [15] . This 
measure quantifies entanglement between the atom and 
different modes in terms of time-dependent distributions. 

The problem of entanglement between an atom and a 
bath of modes, is becoming more interesting when con- 
sidering reservoirs with a spectral gap in their densities of 
states. For such systems, it is well known that an atom- 
photon bound state can be formed [9|, Il6l422j . In view 
of this result, it is reasonable to expect strong quantum 
correlations and entanglement between the atom and the 



reservoir. 

Motivated by this we consider here a two-level atom 
coupled to a model reservoir with a single frequency gap 
in its density of modes. Exploring the dynamics at dif- 
ferent coupling regimes, we are able to show that when a 
trapping state is formed, permanent correlations are ob- 
served. Using the pseudomodes technique [23|, 0[, and 
a tripartite entanglement measure, the tangle [25j, we 
quantify and study the properties of entanglement. Fur- 
thermore, a careful analysis reveals that in the long time 
limit and when a trapping state is formed, the reservoir 
spectrum is not constant in time. This is due to a contin- 
uous coupling between the atom and individual modes, 
which has zero net energy flow, but induces a permanent 
effective coupling between the reservoir modes. In terms 
of the pseudomode description, the population trapping 
arises because of the dark state between the atom and 
one of the pseudomodes. 

This paper is organised as follows. In section [III we 
introduce the model and the pseudomodes method. In 
section HIH we discuss the formation of the trapping state 
and the reservoir dynamics in the long time limit. In sec- 
tion |IV1 an analysis of entanglement dynamics in terms of 
the tangle and the density of entanglement is presented. 
We conclude in section [Vj and in appendix [X] a synopsis 
of the pseudomodes method is provided. 



II. MODEL 



The system we consider in this work, consists of a two- 
level atom coupled to a reservoir of harmonic oscillators 
with annihilation and creation operators a\ and re- 
spectively. Within the rotating wave approximation the 
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Hamiltonian reads (h = 1) 

A 

+ (4l°a)(la|+a A |l a )(0 a |) 



BAND-GAP RESERVOIR 



(1) 



where g\ is the coupling between the mode A and the 
atomic transition |1 ) — |0 o ). The atomic transition 
frequency is ujo whereas the A- mode frequency is uj\. 

For the purposes of the analysis that follows, it is very 
useful to introduce the reservoir structure function D(ou\) 
which reflects the properties of the density of modes [23| . 
This is defined through 



p x {g x f = gD( WA ), 



and is normalized such that 



/oo 
dujD(uj) = 2tt. 
-co 



(2) 



(3) 



With this normalization a measure of the overall coupling 
strength is which is given by 



n 2 



(4) 



In Eq. (|2]) p\ is the density of modes i.e. the number of 
modes with frequencies in the interval uj\ to uo\ + duj\. 

Previous studies revealed that the formation of an 
atom-photon bound state is plausible, when an atom is 
coupled to a reservoir with a gap in its structure function 
[l, [3-22]. It has also been suggested that the formation 
of such a bound state is an indication of entanglement 
between the atom and its environment [26|, [27j . In order 
to explore entanglement dynamics between an atom and 
a reservoir with a gap at a given frequency co> c , we utilise 
the following structure function for the reservoir 



D(lj) = Wt- 



-W 2 



(u> - uj c f + (rx/2) 2 ~'(uj- uj c f + (r 2 /2) 2 ' 

(5) 

This superposition of Lorentzians with the same centre 
frequency co> c , widths Tj, amplitudes Wj and opposite 
signs will result in a gap, i.e. D(co> c ) = 0, if TiW^ = 
^Wi- Because of the normalisation condition (|3j) we 
also have that Wi - W2 = 1. 

Starting with the atom initially excited and the reser- 
voir in the vacuum state, one has to solve the Schrodinger 
equation to obtain the system dynamics for t > 0. 
This can be done either with analytical methods, e.g. 
the Laplace transform [2(J [HI , or numerical integration 
[il, [HI, [l7| . An alternative app roach is that offered by the 
pseudomodes method [23|, [24{ . 

According to this method the reservoir modes are re- 
placed by two degenerate pseudomodes [23|,[23, see figure 
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FIG. 1. Diagrammatic representation of the atom- 
pseudomodes system. The band-gap reservoir is represented 
by two interacting pseudomodes PMi and PM2. The two 
pseudomodes are also coupled to two independent Markovian 
reservoirs, that induce the decay of the pseudomodes at rates 
T[ and r' 2 respectively. The atom couples only to one of the 
two pseudomodes. For a perfect gap i.e. D(uj c ) = 0, the decay 
rate for the first pseudomode T[ is zero. 



[TJ The two pseudomodes are interacting with each other, 
while one of them is also coupled to the atom. Finally the 
two pseudomodes decay at rates T[ and Y' 2 respectively. 

The dynamics of the system are described by a Marko- 
vian master equation [23|, [24[ 



r' 



pit) = - i [Ho, p(t)] - E y ( a jM*) 



(6) 



2a n 



jp{t)a)+ p(t)a)a^ 
with the Hamiltonian 

Hq =uj \l a )(l a \ +o; c (a\ai + a\a^ 

+ fi (4|0a><la|+a 2 |la><0l|) (7) 

+ V (a[a2 + aia^j , 

where a\ (a\) and a 2 (aj) are the annihilation (cre- 
ation) operators for the two pseudomodes respectively. 
The vacuum and excited states for the atom are |0 a ) 
and |l a ). The coupling f^o is given by Eq. (j4j), and 
V = vWiW 2 (ri - T 2 )/2. The two decay rates are 
ri = W 1 T 2 - W 2 T 1 and Y' 2 = - W 2 T 2 . For a 

perfect gap D(uj c ) = 0, the decay rate for the first pseu- 
domode is r; = 0, Y' 2 = (T 1 + T 2 ) and V = vTITJ/2. 
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The solution for the master equation © reads 

P (t) = n i (t)|o a Oi0 2 )(o a Oi0 2 | + m)){i>(t)\, (8) 

where 

|-0(*)> = c a (t)|l a 0i0 2 ) + ai(t)|0„li0 2 ) + O2(t)|0„0il2>. 

(9) 

The vacuum state |0 a 0i02) population is Ilj(t), and the 
probability amplitudes for the atom, and the first and 
second pseudomodes are c a (i), a\(t) and a,2 (i) respec- 
tively. A synopsis of the pseudomode method is provided 
in the appendix [Aj along with the expressions for the 
probability amplitudes c a (i), ai(t), a2(t) and the popu- 
lation Hj(t). 



III. POPULATION TRAPPING AND 
NON-STATIONARY ENVIRONMENT 

For a resonant system uj c = cjo, and when the per- 
fect gap condition T1W2 = F2W1 is met, we have that 
D(uu c ) = and in the long time limit a trapping state is 
formed [24|. Upon solving the equations for c a (t), a\(t) 
and a,2(t), see appendix [Aj and taking their limits for 
t 00 we have that 



and 



c a (oo) = (1 + 77 2 )" 



ai(oo) 7/(1 + rfy 



(10) 



(ii) 



where rj = 2^o/vTiF2. The probability amplitude for 
the second pseudomode is (22(00) = and the population 
of the vacuum state is 



n j (oo)=r/ 2 (l + r/ 2 )- 



(12) 



A plot of \c a (t)\ 2 , |ai(£)| 2 , \a 2 (t)\ 2 and U^t ) for T x = 
10^o an d F 2 = 0.2^o is shown in figure [2(a)] 

From the above three equations it is evident that in the 
long time limit, and in the weak coupling regime rj <C 1, 
a fraction of the population will remain trapped in the 
excite d ato mic state and the the second pseudomode, see 
figure [2(b) , The remaining population is irreversibly lost 
to the reservoir (or more precisely, t o the Markovian part 
of the reservoir (28)). From figures [2(b)] and Eqs. (fTO]) - 
([T2]) we can see that population trapping, i.e. |c a (oo)| 2 , is 
significant for rj < 1. The population lost to the reservoir, 
i.e. the sum of the populations for the pseudomode 1 
l^i (00) 1 2 and the vacuum state Ilj(t), remains low for 



77 <C 1, see Fig. 2(b) As we move to the strong coupling 
regime 77 > 1 losses increase, and eventually for 77 ^> 1 all 
the population is transferred to the reservoir. 

An interesting feature of the trapping state, is that 
in the long time limit the reservoir modes do not reach a 
steady state. This can be evidenced in the reservoir spec- 
trum for t —> 00. Using the definition for the reservoir 
spectrum [29] 
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FIG. 2. (Color online) (a) The populations |c a (t)| 2 (black), 
|ai(£)| 2 (red short-dashed), |«2 (t) | 2 (blue dashed) and Hj(t) 
(green long-dashed), for Fi — 10^o, T2 — 0.2Qo and Wi = 
50W 2 . (b) The final populations |c a (oo)| 2 (black), |ai(oo)| 2 
(red short-dashed) and IIj(oo) (green-dashed) as functions of 
the dimensionless parameter 77 = 2f2 / VTiIV 



and Eq. (|A17p for t(T 1 + T 2 ) > 1 we get the following 
expression for S((J\,t) 



S{u> x ,t) 



7r(4r 2 + ft 2 ) 2 



2S X 



e Jxt/2 sin 



Sxt 
2 



AO 2 (2T - iS x ) 



4(T-i5 x ) 2 + n 2 



(14) 



S(u x ,t) =px\cx(t)f 



(13) 



where the width T and the Rabi frequency Q are given 
in Eqs. (jA14j) and (jA15j) . and S x = oj x - oj c . 

Thus in the long time limit, although the total ex- 
citation in the reservoir is constant, the modes remain 
coupled to each other. As a result of this the population 
distribution between the modes changes, see Fig. [3l The 
oscillatory exchange of population between the modes 
is more pronounced in the weak coupling regime ^0 <^ 
VTiT2, Fig. 3(a)| and is negligible for the strong cou- 
pling regime, Fig. |3(b)| Snapshots of the reservoir spec- 
trum for times t(Ti + T 2 ) ^> 1 are shown in Figs. |3(c) 
and |3(d)| It is also interesting to note that in the weak 
coupling regime displaying the population trapping, the 
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frequency gap imposes strong oscillations in the mode 
populations compared to the single Lorentzian structure 
function case[15], see Fig. 3(c)| In contrast with strong 
coupling and no population trapping, there is a strong 
resemblance in the mode popul ations between the gap 
and single Lorentzian cases (Fig. |3(d)| . 

In order to explore further these features, we plot in 
Fig. |4(a)| the probability current between the atom and 
the A mode KM Ell 



J Xt a(t) = 2lm{p x g x cl(t)c a (t)e^ t }, 
and in Fig. |4(b)| the net probability current 

/oo 
doj x J^ a (t). 
-oo 



(15) 



(16) 



The long time limit, i.e. when t — » oo, is of particular 
interest. We see that while the net probability current 
Fig. |4(b)| approaches zero, in the long time limit, the 
individual reservoir frequency components seen in Fig. 

a) do not decay, but continue oscillating. This does 
not happen in the case of a Lorentzian reservoir coupling 
and appears to be a feature of population trapping in a 
photonic band-gap structure. Subsequently, there exists 
an effective, atom-mediated, coupling between the modes 
even though the atom has reached a steady state. 

For reservoirs with a single Lorentzian structure func- 
tion, the pseudomodes method has provided an intuitive 
insight into memory effects [28|. When an atom is cou- 
pled to such an environment, slowly decaying oscillations 
between the atom and the reservoir are observed in the 
strong coupling limit [TBI, l28j|. These correlation effects, 
are attributed to a memory part of the reservoir that 
is represented by a single pseudomode. The rest of the 
reservoir acts as a Markovian environment that induces 
a slow exponential decay for the memory part. 

In the current system with a frequency gap in the envi- 
ronment, we have two pseudomodes which store informa- 
tion about the state of the atom. The first one ai(t), i.e. 
the one that forms the trapping state with the atom, is re- 
sponsible for the permanent storage of information. The 
second one d2(t) is responsible for short term storage, but 
eventually gets depleted due to its coupling to the rest 
of the reservoir. Though the atom is directly coupled to 
the first pseudomode only, the interaction between the 
pseudomodes gives rise to the trapping of population by 
forming a dark state for the atom - pseudomode one sub- 
system. 

In general, the population trapping signifies the forma- 
tion of an atom-photon bound state. In view of the strong 
permanent correlation effects that dictate the formation 
of such a state, it is reasonable to expect entanglement to 
be also present. In the following section we use both the 
pseudomodes method, and the recently proposed density 
of entanglement [15], to explore entanglement between 
the atom and the reservoir defined by Eq. (|5j). 



IV. ENTANGLEMENT DYNAMICS 

Identifying and measuring entanglement in multi- 
partite systems presents various complications. Apart 
from the case of a two-qubit system, where entangle- 
ment can be identified both for a pure and a mixed state 
[32[ [33|, multi-qubit entanglement is an open problem 
and to date several measures of entanglement have been 
proposed [H, i, [H, HI S-^. For the analysis that fol- 
lows, we will be using two different measures [lEl |25|. 

The first one, called tangle [25j, is a measure of gen- 
uine tripartite entanglement between three qubits. This 
will be used to explore entanglement dynamics in the 
pseudomodes framework. The second one is the recently 
proposed density of entanglement [15]. This measure is 
appropriate for studying entanglement between an atom 
and the continuum of the reservoir modes. It provides 
valuable information regarding entanglement distribu- 
tion between the atom and the modes and between indi- 
vidual modes. 



A. Tangle 

We start our analysis from Eqs. (|5J) and (|9j), i.e. the 
density matrix for the atom-pseudomodes system. For 
this mixed state, the two pseudomodes can be collectively 
described in terms of a single qubit. The two states for 
this collective qubit are 



\o ps ) = |0i0 2 ), 



(17) 



and 

|i P -> 



(Oi(t)|li0 2 >+O2(t)|0il 2 ». 

(18) 



VM*)I 2 + M*)I 2 

Using these expressions, the state \ip(t)) reads 

|#)) = Ca(t)|la<W + V|aiW| 2 + |a2(t)| 2 |0 o lp S ), (19) 

and the density matrix p(t) becomes 



p(t) = n i (t)|0 a ps )<0 o pi | + |#))(#) | . 



(20) 



We should note here, that the states \0 ps ) and |l ps ) are 
both eigenstates with zero eigenvalues for the reduced 
density matrix pi2(t) = tr a {p} for the two pseudomodes, 
where the tracing is over the atomic states |0 a ) and |l a ). 

Entanglement for this "two-qubit" mixed state can be 
quantified in terms of the concurrence [32|, [33| . This can 
be associated with tangle, a measure of tripartite entan- 
glement for a system of three qubits A, B, and C [25j |. 
The tangle tabc expressed in terms of pairwise concur- 
rences reads 



tabc 



r<2 ^2 ri 



AC 



(21) 



where Cab and Cac are the pairwise concurrences for 
the qubit A with B and C respectively, whereas Ca(bc) 



5 




_ L5 -1.0 _ .5 0.0 o.5 1-0 1.5 -1.5 -1.0 -0.5 0.0 .5 1-0 1.5 

(qj x -qj c )/£L () (oj k -oj c )/CL () 



FIG. 3. (Color online) The reservoir spectrum S(uj\,t) as a function of time, (a) for I\ = 10^o, T 2 = 0.2Q and W\ = 50VK2, 
and (b) for T± = 0.5^o, T2 = O.Ol^o and Wi = 50W2- Figures (c) and (d) are snapshots for the reservoir spectrum for Qot = 50 
and for the parameters of figures (a) and (b) respectively. The red dashed lines in figures (c) and (d) is the spectrum for a 
reservoir with a Lorentzian structure function with Y2 = W2 = 0, Wi = 1 and Fi = 10Qo and Ti = 0.5£lo respectively. 



is the concurrence for qubit A and a qubit (BC) that 
collectively describes qubits B and C. 

The above equation can also be written as an inequality 
i.e. 



Ca(bc) - Cab 



C 2 



AC' 



(22) 



The meaning of these two equations is that entanglement 
between the qubit A and the other two qubits, B and C, is 
manifested through direct entanglement with each qubit, 
thus the two concurrences Cab and Cac , and through a 
three-way (tripartite) entanglement i.e. tabc- 

From the three qubit density matrix Eq. (jSJ, we de- 
rive the reduced density matrices for the atom with each 
individual pseudomode i.e. p a ^\ and p„ 2. Using the con- 
currence for a two-qubit system [H, [33j we obtain the fol- 
lowing two expressions for the concurrence for the atom 
with each pseudomode 



and 



C<l 1 (t) = 4\c a (t)f\a 1 (t)f 



C 2 M) = 4\c a (t)\ 2 \a 2 (t)\ 2 . 



(23) 



(24) 



The final step is to calculate the concurrence for the den- 
sity matrix ([20]) . 



This is the concurrence for the atom and the qubit 
that collectively describes the two pseudomodes. The 
calculation is simple and the concurrence C 2 ^ (t) is 



^a,(12) W 



clAt) 



(25) 



i.e. the tangle for the atom and the two pseudomodes is 
zero. From the definition of the tangle, Eqs. (|2T|) and 
([22]) , and the above result, we conclude that entangle- 
ment between the atom and the pseudomodes is mani- 
fested only through two-way entanglement channels. A 
three way entanglement is completely absent. 

In figures 5(a)| and |5(b)| we plot the concurrences 
as functions of time, for the weak and strong coupling 
regimes respectively. In the weak coupling regime, where 
a trapping state is formed, entanglement at a very early 
stage builds up only between the atom and pseudomode 
two which is responsible for the short term storage of 
information. After reaching a peak, it starts decaying 
where at the same time entanglement between the atom 
and the pseudomode one, which is responsible for the 
long-term storage of information, slowly builds up and 
reaches a steady state. Thus the trapping state is also 
an entangled state between the atom and the reservoir. 
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FIG. 4. (Color online) The probability current J\, a (t) Eq. 
(H4|) (a), for Ti = 10O , T 2 = 0.2Q and Wi = 50W 2 , and the 
total probability current Q(t) Eq. (|15|) for the same parame- 
ters (b). 



FIG. 5. (Color online) The concurrences C a ,i(t) (black solid) 
C a ,2(t) (red short-dashed) and C a> (i2)(t) (blue dashed), for 
Ti '= lOHo, E 2 = 0.2Q (a), and for Ti = 0.5ft , T 2 = 0.01Q 
(b). For both figures Wi = 50W 2 . 



For the strong coupling regime, pseudomode one makes 
a negligible contribution in the entanglement dynam- 
ics. Pseudomode two has a strong contribution for short 
times, where the concurrence C a ^{t) quickly increases, 
and then follows a slowly decaying oscillation pattern. 
These oscillations are the signature of a Rabi splitting 
observed in the strong coup ling r egime [15], see also the 
reservoir spectrum in figure |3(b)| 



B. Density of entanglement 

In order to gain further insight into entanglement dy- 
namics, we need to consider entanglement between the 
atom and each of the reservoir modes. For quantifying 
the distribution of entanglement between the atom and 
the individual reservoir modes, and among the reservoir 
modes, we use the density of entanglement [15j. The den- 
sity of entanglement between the atom and modes with 
frequencies in an interval uj\ to uo\ + dco>A is 



£ A (u\,t) =4\c a (t)\ 2 S(uj x ,t), 



(26) 



and the density of entanglement among the reservoir 
modes reads 



£ R {ujx,uj^t) = 2S(u\,t)S(up,i), 



(27) 



where S(uj\,t) is the reservoir spectrum. 

In terms of these two distributions the total entangle- 
ment or concurrence C 2 (t) for the atom and the reservoir 
modes is defined as the sum of an atom-modes contribu- 
tion 



/oc 
duj x £ A (uj x ,t), 
-CO 

and a reservoir contribution 



(28) 



/CO POO 
dw x / dw„£fl(wA,w M >*)- (29) 
-CO J — CO 



The total entanglement reads 

c 2 (t) = c 2 A (t) + cUt). 



(30) 



For the strong coupling regime, as shown earlier, the 
dynamics are similar to those for an atom coupled to 
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FIG. 6. (Color online) The atom-modes density of entanglement £a(cja,£) as a function of time and the mode frequency oj\ 
(a), and the density of entanglement £r(u>\, cj^, t) between a mode uj^ — uj c + O.lQo and the rest of the reservoir modes as a 
function of time (b). Figures (c) and (d) are snapshots for the mode-mode density of entanglement SR(uj\,uj^,t) for Qot — 10 
and Q t = 30 respectively. For all figures Ti = lOHo, T 2 = 0.2Q and Wi = 50W 2 . 



a reservoir with a Lorentzian structure function. Thus, 
entanglement dynamics will be similar for the atom- 
reservoir system under consideration in Ref. [l5j|. The 
main feature for the reservoir density of entanglement is 
a pronounced Rabi splitting. Furthermore, this splitting 
results in decaying oscillations in the atom-modes density 
of entanglement. 

On the other hand, for the weak coupling regime dy- 
namics are different. The formation of the trapping state 
is associated with a continuous population exchange be- 
tween the atom and individual modes whilst the net flow 
of probability is equal to zero as discussed before. As a 
consequence, both ent anglem ent distributions c hange in 
time, as can be seen in |6(a)| for Ea (j^x > t) and in |6(b)| for 
£r(u\, Co^, t). Both distributions do not reach a steady 
state in the long time limi t. Th is fe ature for £r{(jJx^^) 
is also evidenced in Figs. |6(c) and |6(d)| where the den- 
sity of entanglement for the reservoir modes is plotted for 
different times. 

In contrast to this, due to population conservation and 
the fact that the net population exchange Q(t) in the long 
time limit is zero, the total entanglement between the 
atom and the reservoir C\(t) is constant for t —> oo. The 



same is true for the total entanglement for the reservoir 
modes C\(i) and the total entanglement C 2 (t). In figure 
we plot C\(t), Cft(t) and the total concurrence C 2 (t) 
for the weak coupling regime. 

From this we see that at very early times, a rapid 
build up of entanglement takes place between the atom 
and the reservoir. Entanglement between the reservoir 
modes evolves at a much slower rate. Upon reaching a 
maximum, atom-reservoir entanglement follows a decay 
reaching a steady state at about the same time as the 
entanglement between the reservoir modes does. This 
point in time corresponds to the formation of the final 
trapping state between the atom and the reservoir. 



V. CONCLUSIONS 

In this work we have studied correlations and entan- 
glement for an atom-photon bound state. Such states, 
can be formed when an atom couples to a reservoir with 
a gap in its density of modes. Their main feature is that 
in the long time limit the system reaches a steady state 
where the initial atomic excitation energy, is shared be- 
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FIG. 7. (Color online) The concurrence C\(t) (black solid), 
Cji(t) (red short-dashed) and the total concurrence C 2 (t) 
(blue dashed), for Ti = lOHo, T 2 = 0.2Q and Wi = 50W 2 . 



tween the atom and the reservoir. 

Despite the fact that for such a state no change is ex- 
pected in the long time limit, a careful study of intra- 
reservoir dynamics has revealed that this is not the case. 
For a reservoir with a single frequency gap in its struc- 
ture function, we have shown that in the long time limit, 
the atom exchanges energy with individual modes, and 
though the net energy flow is zero, a permanent effec- 
tive coupling between the modes is induced. As a result 
the reservoir spectrum changes with time and a steady 
energy distribution is never reached. 

The existence of the atom-photon bound state is ex- 
plained, in terms of the pseudomode description, with 
the formation of a dark state by the atom and one of the 
pseudomodes which are both coherently coupled to that 
pseudomode which connects to the rest of the environ- 
ment. In general, the population trapping occurs in the 
weak coupling regime whereas with strong coupling and 
no trapping, also the dynamics within the environment 
begins to resemble the one obtained by single Lorentzian 
distribution. Furthermore, we have studied in the de- 
tail the entanglement dynamics between the atom and 
the pseudomodes, and within the environmental modes, 
high-lighting several qualitative and quantitative differ- 
ences between the weak and strong coupling regimes. 

The results and methods used in this work can be ex- 
tended and applied to systems where the reservoir den- 
sity of modes has a broader gap or a more complicated 
structure. Such evolved reservoir structures can be en- 
countered in photonic crystals, or when considering the 
problem of atom-laser outcoupling from Bose-Einstein 
condensates. 
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Appendix A: Solution of the Schrodinger equation 
and the pseudomodes method 

Starting with an atom initially excited and the reser- 
voir in a vacuum state, the system's wave function for 
t > will be 

^(t)) = Ca(t)\la)\0) + $> A (t)|0 a )|V>A>. (Al) 



The collective vacuum state for all reservoir modes u\ is 

io>=ni°A>. ( A2 ) 



and the state with a single excitation in one of the reser- 
voir modes is 

IV'a) = |1a> II I0fc>- (A3) 

At t = we have that c a (0) = 1 and c\(0) = 0. 

The coefficients c a (t) and c\(t) can be obtained by 
solving the Schrodinger equations 



ic x (t)=g x e^ t c a (t), 



(A4a) 



(A4b) 



where the detuning between the atomic transition and 
the mode A is 5\ = uj\ — ujq. The amplitudes in the 
interaction picture are c a (t) = e 1UJot c a (t) and c\{t) = 
e iuJxt c x (t). 

To derive c a (t) and c\(t) one can numerically integrate 
Eqs. (|A4aP and (|A4b|) using a discretization technique 
0, [H, [I?} or by using the Laplace transform [20|, l22j |. 
A different approach is that offered by the pseudomodes 
method [23l. l24j|. The main feature of this technique, is 
that the infinitely many equations for the reservoir modes 
can be replaced by a finite number of equations. Thus 
the computational effort is substantially reduced. In ad- 
dition to this, the pseudomodes method has provided an 
intuitive insight into non-Markovian dynamics, which are 
observed when an atom strongly couples to its environ- 
ment HH. 

When the reservoir structure function is analytic with 
a finite number of poles in the lower complex plane, Eqs. 
(|A4a|) and (IA4bl) can be replaced by a set of equivalent 
equations [23|, [24J . In this new set of equations the atom 
couples to a finite set of fictitious modes, the pseudo- 
modes, where each of these modes has a one-to-one cor- 
respondence to the poles of D(u). 

For the structure function D(uo) in Eq. ([5]), the anal- 
ysis for arbitrary widths (Ti,r 2 ) and weights (Wi, W2) 
was previously carried out, see Ref. [24]. Here we fo- 
cus only the on perfect gap case i.e. D(uo c ) = 0, where 
the equations for the atomic excitation c a (t) and the two 
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pseudomodes a\{t) and a 2 (t) are [2^ 

\C a (t) =LJoC a (t) + ftoa 2 {t) 



idi(t) =u c ai(t) H ^— -a 2 (t), 



id 2 (t) = ( ^ c - i Fl ^ 2 ) a 2 (t) 



(A5a) 



(A5b) 



(A5c) 



+ fi c a (t) + ^^ai(t). 



These equations can be associated to the following master 
equation 

p(t) = - i [H , p{tj\ - Tl ~^ T2 [a\a 2 p{t) 



(A6) 



2 p(t)a\ + p(t)a£a 2 ), 



- 2a 2 / 
with the Hamiltonian 

H =uo|la)(la| (ajdi +a 2 a 2 ) 



(4|0a)(la|+a 2 |l a )(0 1 |) (A7) 



(al 



a 2 + didj 



where di (d{) and d 2 (d^) are the annihilation (creation) 
operators for the two pseudomodes respectively. 
The solution for the master equation (|A6|) reads 



p(t) = n j (t)|0 a 0i0 2 )(0 a 0i0 2 | + \^(t))(^(t)\, (A8) 
where 

$(t)) = c a (t)|l a 0i0 2 > + ai(t)|0 a li0 2 > + a 2 (t)|0 a 0il 2 >. 

(A9) 

The vacuum state population Hj(t) is given by 

n,(i) = ^^ f*dr\a 2 (T)\ 2 . (A10) 

The Fock states with zero or one excitation for the two 
pseudomodes are |0i) (|0 2 )) and |li) (|1 2 )) respectively. 
From Eq. (|A8|) we see that the atom and the pseudo- 
modes are in a mixed state. 

Equations ()A5a|) - (jA5c|) are linear with time- 
independent coefficients and solutions can be easily 
obtained with the Laplace transform method. With the 
initial population for the atom being c a (0) = 1, and 
both pseudomodes in a vacuum state, ai(0) = a 2 (0) = 0, 
we get for c a (t) 



C a (t) 



rir 2 



2 COS {Y 



(AH) 



and for the pseudomodes 



01 (t) 



iojot 



2 v / rIT^»oe 

(4r 2 + n 2 ) 



2T fnt 
n" sm U 



1 - e~ rt cos 



fit 



(A12) 



and 



02 (t) = 



(A13) 



Here we consider only the resonant case cuo = 00 c . The 
decay rate V and the Rabi frequency Q are 



ri + r 2 



and 



n = - ) /i6fig-(r 1 -r 2 )2. 



(A14) 



(A15) 



Substituting (|A12|) in (|A1Q|) we get for the vacuum 
state population 



n 2 



-2Tt 



4r(4r 2 + ^ 2 ) 4r 



(Ai6) 



rcos(Ot) - § sin(Ot) _ 2F£ 

(4r 2 + o 2 ) 6 



Finally using Eqs. (|A1|) and (|A4b[) we get the amplitudes 
ca(£) for the reservoir modes 



c\(t) 



4r 2 + Q 2 



2<5 A 



e^*/ 2 sin' UA 



2 



4ng(2r-u A ) / 
4(r - i5 x ) 2 + n 2 { 



e is x t-rt cog 



2Q 2 (4(iS x r -T 2 ) + Q 2 ) iSxt _ rt fQt\ 

n(4(r - i<5 A ) 2 + ^ 2 J 

(A17) 

where #a = — is the detuning between the A mode 
and the gap frequency uj c . 
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